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A  Computer  Code  to  Calculate  Emission 
and  Transmission  of  Infrared  Radiation 
Through  Non-  Equilibrium  Atmospheres 


!.  INTRODUCTION 

The  Air  Force  Geophysics  Laboratory  is  in  the  process  of  developing  a  compre¬ 
hensive  computer  code.  NLTE,  for  calculating  the  radiance  due  to  infrared-active 
species  in  the  upper  atmosphere.  In  particular,  the  code  deals  with  atmospheric 
conditions  in  which  local  thermodynamic  equilibrium  (LTE-)  cannot  be  assumed;  that 
is,  conditions  under  which  the  populations  of  vibration- rotation  states  cannot 
(necessarily)  be  given  by  a  Boltzman-distribution  with  the  kinetic  temperature. 
Moreover,  it  does  a  line-by-line  calculation,  rather  than  a  band  calculation,  using 

the  properties  of  individual  rotational  lines  obtained  from  the  AFGL  Atmospheric 

12  3  4 

Absorption  Line  Parameters  Compilation.  ’  ’  '  The  output  is  the  integrated 
radiance  from  each  line  under  consideration.  The  individual  lineshapes  can  also  be 
obtained,  if  necessary. 

At  present,  the  non-LTE  conditions  are  simulated  by  the  use  of  a  vibrational 
temperature  profile.  Work  is  in  progress  to  incorporate  into  the  program  computa¬ 
tions  of  the  direct  physical  mechanisms  that  are  responsible  for  the  populations  of 
the  excited  states.  That  is,  the  excitation  of  infrared-emitting  species  by  sunshine, 
earthshine,  collisions,  and  photochemical  reactions  will  be  used  to  directly  calculate 

(Received  for  publication  6  July  1983) 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  21.) 
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the  populations  of  vibrational  states  rather  than  using  the  vibrational  temperature 
profile  as  input. 

The  looking  geometry  is  assumed  to  be  horizon-to-horizon.  It  is  therefore 
parameterized  by  a  tangent  height,  the  lowest  altitude  along  the  line-of-sight  path. 


2.  FORMULATION 

We  consider  the  path  through  which  a  photon  has  to  pass  to  reach  the  observer 
to  be  made  up  of  a  series  of  slabs.  In  each  of  the  slabs  the  kinetic  temperature,  the 
vibrational  temperature,  and  the  density  of  the  atmospheric  constituents  in  the 
ground  and  excited  states  is  assumed  to  be  constant.  The  altitude  of  the  end  points 
of  each  slab  may  not  differ  by  more  than  1  km  (Figure  1). 


Figure  1.  The  Line-of-Sight  Path  Through  a  Layered 
Atmosphere,  the  enhanced  segment  showing  its 
intersection  with  the  ith  layer  over  a  distance 


At  first  glance,  it  might  appear  that  we  are  constructing  a  model  for  a  one¬ 
dimensional  atmosphere.  This,  however,  is  not  the  case.  One  can  make  the  slabs 
as  small  as  desired— so  small  in  fact  that  the  errors  may  be  determined  principally 
by  the  uncertainties  in  the  input  model  atmosphere.  Nor  is  there  any  requirement 
that  the  two  slabs  corresponding  to  a  given  altitude  on  either  side  of  the  tangent 
height  point  have  the  same  composition. 

In  the  equations  that  we  derive,  we  will  assume  that  lines  for  which  the  intensity 
is  calculated  are  non-overlapping.  That  this  is  a  reasonable  assumption  can  be 
seen  from  the  following  argument.  The  two  closest  lines  in  the  15-4m  band  of  the 
C12C>26  01*0-00^  O  transition  are  Q(21  and  Q(4)  which  are  0.0145  cm  1  apart. 
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Recalling  that  the  Doppler  linewidth  at  200  K,  corresponding  roughly  to  70  km  in 
altitude,  is  about  0.  0005  cm'1,  it  can  be  seen  that  the  separation  between  the  Q(2) 
and  Q(4)  lines  corresponds  to  about  30  linewidths.  Apart  from  the  accidental  over¬ 
lap  of  lines  of  different  bands,  that  may  occur  rather  infrequently,  we  expect  most 
of  the  lines  to  be  spaced  several  linewidths  apart.  The  assumption  of  non-overlapping 
lines  thus  appears  reasonable. 

We  express  the  radiance  in  terms  of  n  ,  the  number  of  photons  (of  frequency 

- 1  v 
v  cm  )  per  wavenumber  per  unit  area  per  steradian  which  are  observed  at  a 

point  z  along  the  line-of-sight  path  per  unit  time.  The  change  in  n^  in  a  path  length 

dz  is  described  by 


dnu 
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(1) 


where  h  and  c  are  Planck's  constant  and  the  speed  of  light,  n£  and  are  the  number 
densities  of  the  lower  and  upper  radiating  states,  B(£  -*u)andB(u  "*  £)  are  Einstein 
coefficients  for  absorption  and  induced  emission,  A  is  the  Einstein  coefficient  for 
spontaneous  emission  and  f  is  the  normalized  lineshape  function  such  that 


QO 


f(v  -  i/Q)  dr  = 


1  . 


(2) 


In  Eq.  (1),  n^  n£  B(£  -*  u)  is  the  term  for  the  absorption  of  incident  photons. 

n^  nu  B(u  -♦£)  is  the  induced  emission  term.  The  factor  (hrQ/c),  which  is  commonly 

absent  in  derivations  of  this  equation,  appears  because  of  the  fact  that  n  is  a  photon 

- 1  ^  - 1 

flux  rather  than  an  energy  flux,  and  because  v  is  in  units  of  cm  instead  of  sec  . 

„  The  Einstein  coefficients5  are  connected  by 

g£  B(£  *♦  u)  =  gu  B(u  -»  £)  (31 


and 


A  =  8 7i  hi/5  B(u  -»  £)  (4) 

where  g£  and  gu  are  the  statistical  weights  of  the  lower  and  upper  radiating  states. 
Through  use  of  Eqs.  (3)  and  (4),  the  definition 


5.  Penner,  S.  S.  (1959)  Quantitative  Molecular  Spectroscopy  and  Gas  Emissivities, 
Addison-Wesley,  Reading,  Massachusetts. 
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and  some  algebraic  manipulation  we  see  that  Eq.  (1)  can  be  cast  in  the  form 
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f(v-vQ) 


(6) 


One  can  readily  see  that  the  (1  -y)  factor  corrects  the  simple  absorption  term 
(which  is  proportional  to  n^  n^)  for  stimulated  emission;  the  other  term  in  the 
brackets  is  the  spontaneous  emission. 

Next  we  connect  the  Einstein  coefficient  to  the  line  strength  which  is  available 
on  the  AFGL  Line  Parameters  Compilation.  *  The  tabulated  line  strengths  S(Tg) 
are  given  for  conditions  of  LTE  at  the  standard  temperature  Tg,  which  is  296  K, 
and  are  related  to  B(f  -»  u)  by 


hv  S(T  > 

B(£  -  U)  =  p  (3t  t  (1  -  exp  <-C2i/0/Ts»  1  (7) 

6.  s 


where  Cg  =  1.4388  K/cm  1  is  the  second  radiation  constant  and  Pj(Tg)  is  the  prob¬ 
ability  of  finding  the  lower  vibration- rotation  state  occupied.  In  general,  P^  =  n^  /n 
where  n  is  the  total  number  density  of  the  specie. 

The  exponential  term  in  Eq.  (7)  takes  into  account  the  stimulated  emission  at 
296  K,  and  in  fact  is  simply  y  evaluated  under  conditions  of  LTE  at  T  . 

We  now  define  the  optical  depth,  T  ,  by 
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(8) 


Since  the  quantities  f,  n^t  and  y  are  (by  assumption)  constant  within  the  confines  of 
each  slab,  the  optical  depth  of  the  ith  slab  is 
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(9) 

(The  subscript  i  is  implied  for  f,  n^,  and  y, )  Expressing  n^  in  terms  of  the  total 

number  density,  this  becomes 


Pj  (T) 


'P7TTT 
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The  term  in  square  brackets  now  represents  the  line  strength  for  conditions  pre* 


vailing  in  the  slab,  including  corrections  for  stimulated  emission,  and  this  times 
f(u  -  Vq)  is  the  conventional  absorption  coefficient  k(^).  *  Pg(T)  must  be  evaluated 
with  the  understanding  that  non-LTE  conditions  prevail,  in  general;  the  same  applies 


to  y.  Both  these  quantities  therefore  actually  depend  on  the  vibrational  temperature 
as  well  as  the  kinetic  temperature  of  the  slab. 

To  complete  the  analysis  one  combines  Eqs.  (6)  and  (8)  to  give 


-  %  f  2  ^0 


This  equation  can  easily  be  integrated  across  a  slab.  The  solution  is 


-At... 


■  Vzi'  « 


,  2  ci/q 
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where  the  slab -dependence  of  y  is  finally  made  explicit  with  the  subscript.  In  Eq. 
(12),  the  first  term  represents  the  absorption  within  the  ith  slab  of  radiance  inci¬ 
dent  upon  it  at  Z.;  and  the  second  term  represents  the  contribution  from  the  ith  slab 
alone.  The  strategy  of  the  program  is  now  clear.  The  boundary  condition  at  the 
far  horizon  is  n^Z.)  =  0  for  all  frequencies.  Using  that  as  the  starting  point,  one 
calculates  the  radiance  at  the  near  end  of  each  slab  from  Eq.  (12)  until  the  observer's 
end  of  the  line-of-sight  path  is  reached.  One  thus  obtains  a  line  profile  for  the 
observed  radiance.  This  profile  is  integrated  over  frequency  to  get  the  radiance 


for  the  line  in  question.  At  the  end  the  result  is  converted  to  units  of 

/ 2  , _ i  _■ 


W/cm  -steradian. 

NLTE  also  calculates  the  integrated  radiance  to  be  expected  if  the  line  is 


thin— that  is,  if  At^  is  always  much  less  than  unity  for  that  line.  The  determination 


of  whether  a  line  actually  is  thin  is  made  by  comparing  this  result  with  the  integral 
ovt  i  the  computed  lineshape.  Agreement  within  5  percent  defines  a  "thin"  line. 

The  program  adds  the  integrated  radiance  from  all  the  lines  it  selects  from  the 
line -file  provided  by  the  user.  It  performs  this  operation  for  each  looking  geometry 
specified,  starting  with  the  lowest  tangent  height.  For  each  tangent  height  after 
the  first,  as  the  computation  proceeds  it  checks  to  see  whether  the  line  under  con¬ 
sideration  has  been  found  to  be  thin  for  paths  characterized  by  a  lower  tangent  height. 
If  it  has,  it  is  also  assumed  to  be  thin  for  the  current  path  and  all  paths  characterized 


V-  /./V 
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by  greater  tangent  heights.  The  lengthy  profile  calculations  are  thereby  avoided 
for  this  line,  and  the  integrated  intensity  is  taken  to  be  the  "thin"  result. 


3.  USE  OF  THE  PROGRAM 

NLTE  has  been  written  in  ANSI-standard  FORTRAN  '77.  It  is  still  in  the  pro¬ 
cess  of  being  developed  and  made  more  efficient,  and  many  changes  are  antici¬ 
pated.  A  listing  appropriate  for  radiance  from  CO2  is  given  in  the  Appendix.  At 
present  NLTE  requires  65,000  words  to  compile  and  execute  on  the  CDC-6600  at 
AFGL.  A  high-optimization  option  should  be  selected  at  compile  time  (OPT  =  2) 
unless  changes  are  being  made  and  the  debug  package  is  required.  Except  for  a 
few  seconds  of  overhead  required  for  initialization,  the  execute  time  is  nearly 
proportional  to  the  number  of  lines  used.  It  is  usually  also  proportional  to  the 
number  of  tangent  heights  selected,  but  if  many  of  the  lines  are  thin  the  execute 
time  can  be  reduced  since  [ewer  computations  are  required.  Selecting  the  Doppler 
lineshape  instead  of  the  Voigt  profile  can  cut  the  run  time  by  a  factor  of  three. 

NLTE  requires  information  from  three  data  files,  one  of  which  is  the 
input  stream.  These  files,  which  are  associated  with  units  1,  2,  and  3,  contain 
miscellaneous  program  input,  an  atmospheric  profile,  and  a  coded  AFGL  line  file, 
respectively. 

Table  1  lists  the  miscellaneous  program  input  on  unit  1.  Three  card-images 
are  required,  and  they  are  read  with  list-directed  reads.  This  means  that  comma 
or  blank  delimiters  separate  the  data  elements,  and  there  is  no  need  to  align  the 
latter  in  particular  columns.  All  variables  of  type  character  must  be  enclosed  in 
quotes  (').  Default  values  result  from  blank  fields.  These  are  defined  by  consecu¬ 
tive  commas  or,  in  the  case  of  the  last  fields  on  a  card,  by  a  single  slash  (/). 
Example  1  gives  a  possible  data  record  for  Unit  1. 

The  first  of  these  three  cards  gives  the  specie  of  interest,  and  three  energies 
which  are  required.  The  latter  are  the  energy  of  the  vibrational  transition  in  ques¬ 
tion,  which  is  needed  to  calculate  y;  the  vibrational  energy  of  the  lower  state, 
needed  for  Pj(T);  and  the  energy  of  the  lowest-energy  vibrational  transition,  needed 
to  calculate  the  vibrational  partition  function.  The  program  does  not  recognize 
defaults  for  any  of  these  quantities. 

The  second  card  identifies  the  lines  which  are  used  to  calculate  the  radiance, 
assuming  they  can  be  found  on  the  line  file.  That  is,  the  required  frequency  range, 
band,  branch,  and  rotational  line  designations  are  specified,  and  all  lines  satisfying 
these  requirements  are  used.  The  frequency  range  is  given  by  VMIN  and  VMAX. 

The  band  is  specified  by  the  upper-  and  lower-state  quantum  numbers,  UST  and  LST; 
the  branch,  by  BR;  and  the  rotational  line  by  NRQL.  Acceptable  values  of  BR  are 
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'P'.  'Q',  'R',  and  'A',  where  'A'  implies  all  branches.  Acceptable  values  of  NRQL 
are  integer  numbers  identifying  single  lines  of  a  branch;  and  999,  implying  all  lines 
of  the  selected  branch(es). 

NDP  is  a  code  specifying  the  lineshape  to  be  used  in  the  computation.  A  value 
of  0  causes  the  Voigt  profile  to  be  used;  1  gives  the  Doppler  lineshape. 

The  program  recognizes  defaults  of  0,  0,  and  50000  for  NDP,  VMIN,  and 
VMAX  respectively. 

The  third  card  gives  the  tangent  heights  which  define  the  looking  geometry  for 
each  desired  case.  TANIN  is  the  lowest  tangent  height,  TANF  the  highest,  and 
INTRVL  is  the  spacing. 


Table  1.  Miscellaneous  Program  Information  on  Unit  1 


Quantity 

Description 

Units 

Type 

Example 

(CARD  1A) 

MOL 

Molecular  formula 

Character 

'C02' 

ISO 

Isotope  code 

cm  "  J 
cm  "  j 
cm  “ 1 

Integer 

626 

VIBE 

Vib  energy  of  tHe  transition 

Real 

667.  3 

VIBL 

Vib  energy  of  lower  state 

Real 

0 

VI BQ 

Vib  energy  for  partition  fn 

Real 

667.  3 

(CARD  IB) 

BR 

Branch  (P,  Q,  or  R) 

Character 

'H ' 

NRQL 

Rotational  quantum  line 

Integer 

999 

UST 

Upper  vib  state  (AFGL  not  'n) 

Chare  cter 

'01 101' 

LST 

Lower  vib  state  (AFGL  not  ’n) 

Character 

'00001' 

NDP 

Lineshape  option 

cm  ’  * 

Integer 

0 

VMIN 

Lowest  wavenumber  desired 

Real,  int 

600 

VMAX 

Highest  wavenumber  desired 

cm  "  * 

Real,  int 

7  50 

(CARD  1C) 

TANIN 

Initial  tangent  height 

km 

Integer 

70 

TANF 

Final  tangent  height 

km 

Integer 

80 

INTRVL 

Examination  height 

km 

Integer 

5 

Example  1  Example  of  an  Input  Data  Recorc 

'C02\  626,  667.3,  0,  667.3 

'R'.  999,  '01101',  '00001'/ 

70,  80,  5 

on  Unit 

1 

Unit  2  contains  the  atmospheric  profile  which  is  to  be  used  in  the  radiance  cal¬ 
culation.  The  file  providing  this  information  is  headed  by  a  single  card-image  con¬ 
taining  identifying  comments.  Each  succeeding  card-image  contains  the  information 


listed  in  Table  2,  in  the  format  given  there.  No  more  than  250  cards,  which 
correspond  to  successively  higher  altitudes  with  a  1-km  spacing,  may  be  read  in. 
The  lowest  altitude  must  be  less  than  or  equal  to  TANIN.  The  highest  must  be 
greater  than  or  equal  to  TANF  +  50,  since  the  calculations  are  done  using  50  layers. 

A  procedure  for  obtaining  the  atmospheric-profile  data  files  in  the  form  re¬ 
quired  by  NLTE  for  input  on  Unit  2  is  described  in  the  next  section 

Unit  3  contains  the  coded  line  file.  This  is  needed  to  determine  the  absorption 
and  reradiation  of  energy  by  each  line  in  the  wavelength  range  under  consideration. 
The  contents  of  such  files,  and  the  means  of  obtaining  them  from  the  complete  data 
base  which  has  been  developed  at  AFGL,  are  given  in  Section  5. 


Table  2.  Atmospheric  Profile  Information  on  Unit  2 


Quantity 

Description 

Units 

Format 

A  (integer) 

Altitude 

km 

15 

TRTMP 

Kinetic  temperature 

Kelvin 

F 10.  3 

ALCOR 

Atmospheric  pressure 

atmos 

E  12.  5 

RHO 

Density  of  radiating  molecule 

cm  *3 

E  12.  5 

VBTMP 

Vibrational  temperature 

Kelvin 

F 10.  3 

4.  A  PROCEDURE  FOR  THE  ATMOSPHERIC  PROFILES 

A  procedure,  called  ATMOS,  is  available  at  AFGL  for  the  purpose  of  allowing 
users  to  create  databases  in  the  format  required  to  input  to  NLTE  on  Unit  2.  The 
data  presently  built  into  ATMOS  are  profiles  very  much  like  those  in  the  1976  U.  S. 
Standard  Atmosphere.  The  vibrational  profiles  can  either  be  read  as  input  to  this 
procedure,  or  derived  in  a  very  simple  (but  contrived)  way  from  the  kinetic  tempera 
ture.  The  alternative  to  reading  vibrational  temperatures  directly  is  to  add  or 
subtract  a  fixed  number  of  degrees  to  or  from  the  kinetic  temperature.  Interpola¬ 
tion  procedures  are  built  in  so  that  input  data  for  the  vibrational  temperature  are 
not  needed  at  intervals  as  closely  spaced  (1  km)  as  they  are  on  output. 

The  procedure,  ATMOS,  is  in  a  very  preliminary  stage  of  development.  In  the 
future,  its  data  base  will  be  expanded  to  include  various  standard  atmospheric  pro¬ 
files  and  some  other  features  as  well.  However,  it  can  be  used  in  its  present  form 
on  the  AFGL  central  computer  system,  as  shown  in  Example  2. 

ATMOS  requires  user  input  at  two  places:  the  BEGIN  card  and  the  data  card 
following  the  end-of-record  in  the  input  stream.  Also,  if  vibrational  profiles  are 
to  be  read  in,  a  logical  file  called  TAPE4  must  be  attached.  The  BEGIN  card  has 
two  parameters  which,  if  selected,  cause  the  output  data  file  to  be  cataloged  as  a 


permanent  file.  The  parameters  are  PFN  and  NAME,  and  are  used  in  the  example. 
It  is  acceptable  to  omit  them  completely,  in  which  case  the  form  is  simply  BEGIN, 
ATMOS,  P,  and  the  data  are  not  saved  on  a  permanent  file.  In  either  case,  the  first 
ATTACH  card  should  be  used  exactly  as  in  the  example. 

The  data  card  should  contain  six  quantities,  listed  in  Table  3.  They  are  read 
with  a  list -directed  read,  discussed  above.  Most  of  them  are  variables  of  type 
character,  and  must  therefore  be  enclosed  in  quotes.  MOL  specifies  the  radiating 
molecule.  Acceptable  values  of  DORN  are  'DAY'  or  'NIGHT',  reflecting  the  fact 
that  some  species'  density  profiles  have  diurnal  variations.  HMIN  and  HMAX  are 
the  minimum  and  maximum  altitudes  of  interest. 


Table  3.  Input  Data  for  the  Procedure,  ATMOS 


Variable 

Description 

Type 

Default 

MOL 

Molecular  Formula 

Character 

DORN 

Day  or  night  profile 

Character 

'NIGHT' 

HMIN 

Lowest  altitude  (km) 

Integer 

0 

HMAX 

Highest  altitude  (km) 

Integer 

250 

CODE 

Code  for  vib  temp 

Character 

'ADD' 

FMT 

Format  for  read,  or 

Character 

'(I5.F10.  3)’ 

number  of  degrees  to  add 

'O' 

EXAMPLE  2  Use  of  the  Procedure,  ATMOS 


USER1.T30.  1234  USER 

ATTACHiP,  PROCFILE,  ID=  WINTERS,  MR=1) 

ATTACH(TAPE4,  USERSFILE,  ID=USER,  MR=1) 

BEGINIATMOS.  P,  PFN=OUTPUTFILE,  NAME=USER) 

EOR 

'C02',  ’NIGHT’.  70,  190,  ’READ’, '(15,  34X,  F10.  3)’ 


CODE  and  FMT  have  different  forms  depending  on  whether  the  vibrational 
temperature  profile  is  to  be  read  from  TAPE4,  or  the  vibrational  temperature  pro¬ 
file  is  to  be  determined  from  the  kinetic  temperature  profile.  Meaningful  values 
for  CODE  are  'READ'  and  'ADD'.  Anything  else  causes  the  vibrational  temperature 
to  be  set  equal  to  the  kinetic  temperature  everywhere. 

If  CODE  is  read  as  'READ',  the  vibrational  temperature  is  read  in  a  user- 
selected  format,  given  by  FMT.  Each  card-image  on  TAPE4,  except  for  a  header 
card  containing  identifying  information,  should  contain  an  altitude  and  a  vibrational 
temperature  in  that  order.  FMT  should  then  specify  two  fields,  one  integer  and  one 
real,  for  reading  these  quantities.  The  default  value  of  FMT  is  '(15,  F10.  3)',  but 
inside  the  quotes  and  parentheses,  which  are  required,  any  format  is  acceptable. 
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If  CODE  is  read  as  'ADD',  FMT  must  be  a  character  representation  of  the 
number  of  degrees  to  add  to  the  kinetic  temperature.  Possible  values  are  '+10', 

'50',  'O',  and  '-50'.  Default  is  'O’. 

5.  THE  AFGL  ATMOSPHERIC  ABSORPTION  LINE  PARAMETERS  COMPILATION 

1  or  many  years,  AFGL  has  maintained  and  continually  updated  its  own  com¬ 
pilation  of  atmospheric  absorption  line  parameters  for  transitions  from  the 

1-4 

millimeter  region  through  the  visible.  These  are  described  in  various  publications. 

The  compilation  is  the  responsibility  of  Dr.  L.  S.  Rothman,  AFGL/OP1. 

For  users  of  the  AFGL  computer  system,  the  data  reside  on  a  disk  pack  in  the 
CDC-6600  computer  center.  For  workers  at  other  installations,  the  compilation 
is  available  on  two  magnetic  tapes.  (It  is  also  accessible  on  two  CC  tapes  at  AFGL, 
but  use  of  the  disk  pack  is  preferable.  )  The  compilation  was  originally  developed 
in  two  parts.  The  "Main"  compilation  contained  data  from  the  principal  sources  of 
atmospheric  infrared  absorption.  These  molecules,  corresponding  to  a  gas  code 
of  1  through  7,  are  HgO,  COg,  O^,  NgO,  CO,  CHj,  and  C>2  respectively.  The  trace- 
gas  file  contained  data  for  the  molecules  whose  codes  are  8  through  28.  ^  In  fact, 
this  division  still  exists  on  the  magnetic  tape  compilations.  A  merged  line  file  was 
recently  written  to  a  single  disk  pack,  however,  making  all  the  data  available  (at 
AFGL)  by  use  of  simple  procedures. 

The  data  exist  in  card-images  on  the  pack  and  tapes,  with  one  card-image 
corresponding  to  each  transition.  There  are  almost  300,  000  transitions  (lines)  in 
all,  arranged  in  order  of  increasing  wavenumber.  They  are  grouped  in  files  spanning 
100  wavenumbers,  each  one  divided  into  records  containing  up  to  250  lines  (pack)  or 
40  lines  (tapes).  On  the  pack  the  files  are  in  buffered  binary  form;  on  the  tapes, 
buffered  BCD. 

The  data  available  for  each  transition  are  given  in  Table  4.  ^  There  are  several 
programs  available  for  accessing  subsets  of  the  complete  compilation  and  making 
them  available  in  particular  formats  for  specific  purposes  (for  example,  input  to  the 
LTE  radiance  program,  FASCODE  ).  The  formatting  listed  in  Table  4  is  the  output 
specification  of  a  procedure,  WRITE,  written  by  Dr.  Rothman  to  obtain  a  coded  line 
file  from  the  disk  pack.  Generally,  subsets  of  the  complete  data  set  are  located 
according  to  particular  (1)  frequency  ranges,  (2)  molecules,  (3)  isotopes,  and 

6.  Clough,  S.  A.  ,  Kneizys,  F.X.  ,  Rothman,  L.S.,  and  Gallery,  W.  O.  (1981) 

Atmospheric  Spectral  Transmittance  and  Radiance:  FASCOD1B, 

AFGL-TR-8 1-0269,  Ad  A104832;  also  Sf’IEl  2^7  (Atmos  Transmission) 

152-166. 
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(4)  vibrational  bands.  However,  by  reading  the  line  file  and  checking  any  of  the 
parameters  given,  one  can  make  a  selection  based  on  other  criteria  as  well. 


Table  4.  Line-File  Data  for  Each  Transition 


Format 

Symbol 

Line-File  Datum 

F 10.  4 

V 

Resonant  frequency  (cm-ll 

E10.  3 

s 

Line  intensity  at  296  K  (cm- 1 /(molecule-cm2)) 

F5.4 

a 

Halfwidth  (air  broadening)  (cm-l/atmos) 

F 10.  3 

E" 

Energy  of  the  lower  state  (cm-1) 

2A8.A10.A9 

Quantum  numbers,  line  identifiers 

13 

Entry  code  for  these  data 

14 

ISO 

Isotope  code 

13 

MOL 

Gas  code 

The  meaning  of  the  physical  data  (frequency,  and  so  on)  listed  in  Table  1  is 
self-evident.  The  quantum  numbers  are  discussed  in  Reference  1  (and,  to  a  small 
extent,  below).  The  entry  code  is  usually  irrelevant  for  the  user.  The  gas  code 
is  an  integer  between  1  and  28,  as  mentioned  above.  The  isotope  code  is  also 
explained  in  Reference  1.  It  is  a  string  of  integers  representing  the  last  digit  in 
the  atomic  weight  of  each  atom,  given  in  the  order  in  which  the  atoms  exist  in  the 
molecule.  For  example,  the  main  isotope  of  CC>2  is  identified  by  the  code  626, 
since  the  molecule  consists  of  atoms  of  atomic  weight  16,  12,  and  16. 

5.1  Obtaining  Files  From  the  Disk  Pack 

The  procedure  WRITE,  mentioned  above,  is  the  easiest  way  for  users  at  AFGL 
to  obtain  data  from  this  compilation.  Example  3  gives  a  control-card  sequence 
which  obtains,  catalogs  on  permanent  file,  and  lists  all  ozone  lines  between  2000  and 
2020  cm  .  For  any  job  using  this  procedure,  the  first  three  cards  must  be  pre¬ 
pared  exactly  as  in  the  example,  except  for  the  user's  name  and  problem  number, 
the  specification  of  the  logical  file  name  (P)  and  possibly  the  time  allocated  to  the 
job.  The  BEGIN  card  contains  optional  parameters,  described  below,  which  deter¬ 
mine  the  disposition  of  the  information  acquired  from  the  compilation.  The  data 
card(s)  following  the  end-of- record  (EOR)  provide  the  information  the  procedure 
requires  to  identify  the  desired  lines. 

The  information  to  be  taken  from  the  data  cards  is  listed  in  Table  5.  Each  card 
is  read  with  a  list-directed  read,  which  means  that  the  data  fields  (eight  in  number) 
should  be  separated  by  comma  delimiters  and  no  attention  need  be  paid  to  the  col¬ 
umns  used  for  each  field.  Default  values  result  from  blank  fields  (consecutive 
commas).  After  the  last  field  containing  a  non-default  value,  a  slash  (!)  can  be  used 
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to  terminate  data  input  for  that  card.  Either  the  slash  or  a  sufficient  number  of 
commas  to  define  each  field  must  be  used.  All  input  data  of  type  character  must 
be  surrounded  by  quotes,  like  '03'  in  Example  1.  Leading  blanks  inside  the  quotes 
will  give  unpredictable  results.  Anywhere  else,  blanks  are  irrelevant. 

EXAMPLE  3  Use  of  the  Procedure  WRITE 

USER1,  T30,  CM50000,  STMFA,  PK.  1234  USER 

ATTACH(P.  PROCEDURES,  ID=  ROTHMAN,  MR=1) 

MNT(VSN=FAS53,  SN=FASPK> 

BEG  IN (W  RITE,  P,  PFN=OZONELINES,  NAME=  USER,  COPFIL=YES) 

EOR 

2000,  2020,  '03', , , , ,  'PFILE' 


Table  5.  Input  Data  for  the  Procedure,  WRITE 


Datum 

Quantity 

Type 

Default 

VI,  V2 

Frequency  range  (cm-1) 

Integer  or  real 

MOL 

Molecular  formula 

Character 

All  molecules 

ISO 

Isotope  code  (Reference  1) 

Integer 

All  isotopes 

UST,  LST 

Band  quantum  numbers 

Character 

All  branches 

SCRIT 

Threshold  strength 

Real 

All  strengths 

DISPOS 

File  disposition  code 

Character 

List  only 

VI  and  V2  must  be  specified,  as  there  is  no  default.  MOL  is  an  alphanumeric 
symbol  giving  the  molecular  formula.  ISO  is  the  isotope  code  explained  earlier,  UST 
and  LST  are  the  upper-  and  lower-state  quantum  numbers  of  the  vibrational  band 
desired.  The  quantum  numbers  are  different  for  different  molecules,  and  their 
meaning  is  discussed  in  Reference  1.  For  the  purposes  of  WRITE,  one  enters  con¬ 
secutive  digits  inside  the  quotes,  for  example  'll  101'. 

SCRIT  is  an  intensity  threshold.  Lines  which  are  weaker  than  this  threshold 
are  ignored.  Because  of  the  small  magnitudes  involved,  an  E  descriptor  should  be 
used. 

The  data  selected  by  the  program  can  either  be  listed  on  OUTPUT  or  written  to 
TAPE8,  which  is  later  made  into  a  permanent  file.  DISPOS  is  a  code  used  to  select 
the  desired  disposition.  Entering  either  'PFILE'  or  'PUNCH'  in  this  field  causes 
TAPE8  to  be  written.  Anything  else  invokes  the  list  option.  One  can  list  some  data 
and  save  other  data  by  entering  different  disposition  codes  on  different  cards. 

Lines  satisfying  the  requirements  of  each  separate  data  card  are  frequency- 
ordered  on  whichever  file  they  are  directed  to,  OUTPUT  or  TAPE8.  However,  lines 
from  one  data  card  will  always  follow  those  from  preceeding  cards.  If  the  frequency 
ranges  are  increasing  and  do  not  overlap,  frequency-ordering  is  preserved. 


Otherwise,  for  TAPES  only,  a  SORT/MERGE  option  (see  below)  can  be  specified  on 
the  BEGIN  card  to  restore  frequency-ordering. 

The  BEGIN  card  invokes  the  procedure  WRITE  (which  is  found  on  logical  file  P 
in  our  example).  It  allows  up  to  four  optional  parameters  which  are  used  for  three 
separate  purposes:  to  preserve  the  frequency-ordering  on  TAPE8,  as  mentioned 
above;  to  give  the  permanent  file  name  and  user  identification  required  for  cataloging 
TAPE8;  and  to  cause  the  information  on  TAPE8  to  be  copied  to  OUTPUT.  The  form 
of  each  of  these  parameters  is  KEYWORD^  VALUE,  where  KEYWORD  is  fixed  and 
VALUE  is  supplied  by  the  user.  Table  6  describes  these  parameters.  PFN  and 
NAME  should  be  either  both  be  specified,  or  neither.  The  parameters  can  be 
entered  in  any  order  following  the  logical  file  name.  Parameters  which  are  omitted 
assume  default  values,  the  effect  of  which  is  given  in  Table  6. 


Table  6.  Parameters  for  the  Procedure,  WRITE 


Keyword 

Value 

Default 

SORTM 

YES 

No  SORT/MERGE 

PFN 

Permanent  File  Name 

No  catalog  is  done  even 
if  PFILE  is  specified 

NAME 

User's  ID 

on  some  data  cards 

COPFIL 

YES 

No  list  of  TAPE8 

It  is  possible  to  invoke  WRITE  more  than  once  in  a  single  job,  with  different 
parameters,  so  long  as  TAPES  is  returned  after  each  usage  and  so  long  as  one  data 
record  is  provided  for  each  usage. 

The  WRITE  output  specification  given  in  Table  4  is  that  used  on  TAPE8.  It  may 
change  slightly  for  certain  frequency  ranges,  since  sometimes  more  significant 
figures  are  available  than  this  particular  format  allows.  The  field  widths  are  always 
the  same,  however.  On  OUTPUT,  when  the  list  option  is  selected,  the  fields  are 
spread  out  for  tase  in  reading  and  to  make  space  for  annotation. 

Example  4  contains  a  small  subset  of  lines  from  the  11101-00001  band  of  COg. 
in  the  format  assigned  by  WRITE  and  described  above.  Example  5  gives  a  possible, 
albeit  unusual,  set  of  requirements  for  linefiles,  with  many  different  forms  used 
on  the  data  cards.  Only  three  sets  of  lines  are  written  to  TAPE8,  and  these  are 
then  ordered  before  being  cataloged  and  copied.  The  second  call  to  WRITE  pro¬ 
duces  the  lines  of  which  Example  4  is  a  subset. 
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EXAMPLE  4  Portion  of  a  Coded  File  from  WRITE 


2075. 2960 

7.  868E-24.  0770 

2.341 

11101 

00001 

p 

2 

482 

626 

2 

2076.8632 

3.  363E-23.  0770 

2.  341 

11101 

00001 

Q 

2 

482 

626 

2 

2076.  8788 

5.  894E-23.  0770 

7.  804 

11101 

00001 

Q 

4 

482 

626 

2 

2076.  9034 

8.  166E-23.  0760 

16. 389 

11101 

00001 

Q 

6 

482 

626 

2 

2076. 9369 

1.  009E-22.  0750 

28.  095 

11101 

00001 

Q 

8 

482 

626 

2 

2076.  9793 

1.  160E-22.  0750 

42.  922 

11101 

00001 

Q 

10 

482 

626 

2 

2077. 0308 

1.  265E-22.  0740 

60.  871 

11101 

00001 

Q 

12 

482 

626 

2 

2077. 0912 

1.  325E-22.  0740 

81.  940 

11101 

00001 

Q 

14 

482 

626 

2 

2077. 1607 

1.340E-22.  0730 

106. 130 

11101 

00001 

Q 

16 

482 

626 

2 

2077. 2393 

1,  316E-22.  0720 

133.439 

1 1 101 

00001 

Q 

18 

482 

626 

2 

2077. 3270 

1.  258E-22.  0720 

163. 868 

11101 

00001 

Q 

20 

482 

626 

2 

2077.4239 

1.  173E-22.  0710 

197. 416 

11101 

00001 

Q 

22 

482 

626 

2 

2077. 5299 

1.  069E-22.  0710 

234. 083 

11101 

00001 

Q 

24 

482 

626 

2 

2077.6373 

1.  252E-23.  0780 

0.  000 

11101 

00001 

R 

0 

482 

626 

2 

EXAMPLE  5  Possible  Usage  of  the  Procedure  WRITE 

USER2,  T30,  CM50000,  STMFA,  PK.  1234  USER 

ATTACH(P,  PROCEDURES,  ID=  ROTHMAN,  MR=1) 

MNT(VSN=FAS53,  SN=FASPK) 

BEGIN(W RITE,  P.  SORTM=YES,  PFN=MISCLINES,  NAME=  USER,  COPFIL=YES) 
RETURN(TAPE8) 

BEGINIWRITE,  P,  PFN=C02LINES,  NAME=USER) 

EOR 

1220.000.  1221.000/ 

1220.  1221, 'H20'/ 

1220  ,  1221,  'H20  '  .181  / 

1220.  1221, ,, 1 100',  '000'  / 

1220,  1221.....5.E-22  / 

1220.  1221, ....  5.  E-22,  'PFILE' 

1220.  1221,, ,  '000001  IP,  00000000',,  'PUNCH' 

1220.,  1221.,  211,,,, 'PFILE'/ 

EOR 

2000,  2173,  'C02',  626,  '11101',  '00001', ,  ’PFILE’ 


5.2  Obtaining  Files  From  the  Magnetic  Tapes 


Users  who  do  not  have  access  to  the  disk  pack  can  obtain  the  line-file  data  from 
the  two  magnetic  tapes.  The  tapes  contain  files  spanning  100  wavenumbers  and  each 
file  is  split  into  records  containing  up  to  40  transitions  (the  last  record  generally 
being  shorter).  Each  card-image  consists  of  eight  60-bit  words,  and  there  is  a 
single  word  at  the  beginning  of  each  record  specifying  the  number  of  transitions  on 
the  record.  The  typical  record  therefore  consists  of  321  words. 

The  data  on  the  tapes  are  the  same  as  is  given  in  Table  4,  and  in  fact  the  format 
appropriate  for  decoding  the  record  is  exactly  that  given  in  Table  4  except  for  the 
35  characters  specifying  the  quantum  numbers  and  line  identifiers.  The  format  for 
this  portion  of  the  card-image  depends  on  the  molecule  in  question  and  is  described 
in  Reference  1.  In  this  context  (magnetic  tapes)  the  ordering  and  formatting  given 


in  Reference  1  is  correct,  unlike  that  of  the  new  disk  pack  where  uniformity  dic¬ 
tated  certain  changes. 

In  order  to  obtain  a  line  file  the  coded  records  containing  the  data  must  some¬ 
how  be  read  into  arrays,  and  the  data  must  be  checked  to  see  if  they  are  included 
in  the  user's  desired  subset.  In  other  words,  a  program  to  serve  the  function  of 
procedure  WRITE  is  required.  On  CDC  systems,  each  record  can  be  read  with 
a  BUFFER  IN  statement.  One  can  then  use  the  internal  file  feature  of  FORTRAN  '77 
as  a  means  of  converting  character  information  to  numerical  information.  This 
feature  utilizes  a  variable  or  array  of  type  character  as  the  internal  file.  For¬ 
matted  read  and  write  operations,  to  and  from  this  region  of  memory,  perform 
the  conversion.  This  approach  is  superior  to  the  use  of  DECODE  for  several 
reasons,  one  of  which  is  that  it  is  ANSI- standard  and  hence  available  to  all 
FORTRAN  '77  users. 

Example  6  gives  a  small  portion  of  a  possible  FORTRAN  program  which  could 
be  used  to  examine  the  line-file  data  on  a  CDC  system.  The  tape  is  identified  with 
TAPE5.  After  the  record  is  buffered  in,  the  first  word,  A(l),  is  checked  to  deter¬ 
mine  the  number  of  transitions  on  the  record.  The  internal  file  used  here  is  the 
variable  C.  Then,  using  array  B  as  the  internal  file,  the  data  for  all  transitions 
are  read  into  arrays  according  to  FORMAT  statement  102.  All  the  information  is 
now  in  numerical  form  and  can  be  manipulated  with  normal  arithmetic  operations, 
except  the  35  characters  of  identifying  information  which  are  left  in  character  form 
and  can  therefore  be  compared  only  with  other  character  variables. 

The  loop  DO  300  performs  simple  comparisons  for  each  line  to  see  if  the  data 
satisfy  the  user's  criteria.  The  third  and  fourth  IF  statements  compare  substrings 
of  the  Ith  element  of  CH  with  the  codes  for  the  desired  branch,  defined  in  the  DATA 
statement. 

This  sketchy  example,  with  its  data  and  formats  "hardwired"  in,  is  obviously 
insufficient  to  duplicate  the  capability  of  procedure  WRITE,  but  the  basic  idea  can 
be  developed  into  something  considerably  more  general.  Note  that  there  is  no  need 
to  dump  all  of  the  35  characters  of  identifying  information  into  a  character  array. 

Some  of  that  information,  such  as  the  upper-  and  lower-state  quantum  numbers,  is 
in  the  form  of  digits  and  could  be  converted  to  integers  for  arithmetic  comparisons. 
Alternatively,  any  or  all  of  it  could  be  ignored  altogether. 

This  also  neglects  the  matter  of  file  positioning — that  is,  skipping  files  and 
records  until  the  record  containing  the  lowest  desired  frequency  is  in  position  to  be 
read.  This  is  obviously  desirable,  to  avoid  reading  and  simply  discarding  hundreds 
or  thousands  of  records.  Since  files  span  100-wavenumber  regions,  the  number  of 
files  to  skip  is  easily  determined.  For  example,  to  get  lines  in  the  region  650-750  cm 
six  complete  files  are  skipped:  SKIPF(TAPE5, 6,  17,  C).  If,  in  addition,  the  lines 
at  650  cm  1  are  known  to  be  in  the  (N+l)st  record,  N  records  can  then  also  be 


skipped  within  the  seventhfile):  SKIPF(TAPE5,  N,  1,  C).  Otherwise,  as  is  generally  the 
case,  one  must  read  each  record  of  the  correct  file  to  discover  its  range  of  wavenumbei 
For  users  at  installations  without  CDC  equipment,  the  procedures  for  reading 
the  coded  tape  are  likely  to  be  quite  different.  In  particular,  the  different  word 
sizes  and  character  codes  have  to  be  accounted  for. 

EXAMPLE  6  Portions  of  a  possible  FORTRAN  code  for  reading  magnetic 
tape  line  files. 

PROGRAM  X 
CHARACTER-80  B(40) 

CHARACTERS  CH(40) 

CHARACTE R*10  C,BRDES(2) 

DIMENSION  A(32 1),  V(40),  STS(40),  ALF(40),  EDP(40),  ISO(40),  MOL(40) 

DATA  BRDES/ '  01  101’, '00001'/ 

OPEN(5) 

1 00  FORMA  T(8A  10) 

101  FORMAT(IIO) 

102  FORMAT(F10.  4,  E10.  3,  F5.  3,  F10.  3,  A35,  3X,  14,13) 

MDES  =  2 
IDES  =626 

150  BUFFER  IN  (5,0)  (A(l),  A(321)) 

IF(UNIT(5))200,  900,  900 
C 

C  CHECK  N.  THE  NUMBER  OF  LINES  ON  THIS  RECORD 
C 

200  WRITE(C.  100)A(1) 

READ(C,  101  )N 
C 

C  DECODE  THE  REST  OF  THE  RECORD 
C 

M  =  8*N  +  1 

WRITE(B,  100)(A(I),  1=2,  M) 

READ(B,  102)(V(I),  STS(I),  ALF(I),  EDP(I),  CH(I),  ISO(I),  MOL(I),  1=1,  N) 

C 

C  CHECK  FOR  THE  DESIRED  MOLECULE  (MDES),  ISOTOPE  (IDES)  AND 
C  BRANCH  (BRDES) 

C 

DO  300  I  =  1,  N 

IF(MOL(I).  NE.  MDES)GO  TO  300 
IF(ISO(I).  NE.  IDES)GO  TO  300 
IF(BRDES(1).  NE.  CH(I)  (3:12))GO  TO  300 
IF(BRDES(2). NE.  CH(I)  (18:27)GO  TO  300 
C 

C  STORE  THE  DATA  FOR  THIS  LINE 


300  CONTINUE 
GO  TO  150 


ERROR  BRANCHES 
900  CONTINUE 


i  \T 
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PROGRAM  KITE 


INTEGER  I  AN  IN « TANF  , A .MCL.Ht S, GEGV. FL AG 


CHARACTER**  XX  (T)  ,  HCL 


I J  1:#  T  * 


no  n*.i  nikfiKff.wtfi.il 


OIMENSICK  K  VZ  (  101  ,  4)  ,  OVV<  1C  1) 

OIREAXIQN  I  HINT  125  il  xtf£  TNP  (25E1  »RHC  (  E5J11 ,  CRAT  125  G1  ,  AlCGR  (  25  Cl 
OINENSICN  IZt2'3,4),G8e  125C1.SI250  I, CHIC  12*0), TGT125C1 
DIMEKS1QK  J.H541il  .?0LL*4iii 

OINENSICN  V  ( 5  4  £ ) ,  S  T  S  ( 5  4  6 )  »ALF  (546)  ,E0»  15  16) 

DIMENSION  5  AO  (41 .SLHRAC  <41  .HI?(4>  .  1HT4BIU  ,F1  An  I  1.1 _ _ 


13 

C 

lH 

_  OF£N(l,FUc  =  *IKPUT ') 

15 

0PEN(2) 

_14 

_  OPEN! 3) 

1Z 

OPEN{4,FIL£x»ClTPUT‘ 

1  FORMAT  <//) 
3_  FORMAT  (IX) 


_Q££IN£_ALL  -C0M1TAN  ISLAJ'D— 1NIT IAL  VALLES 
EARTHS  x  EARTH'S  RACIUS  (KM 

IS  =  S  T  A  NO  A  R  0  LINE  STREKGTI-  IEMFERATLRE  ( K )  FCS  AFGL  TAPE 
BOITZ  =  ROLIZNAN'S  CONSTANT  (cPGS/DEG) 

8K  =  eClTZMAN'S  FACTOR  (CM-l/DEG).  THIS  FACTOR  IE  NEEEEC  IN 
EXFRESSIONE  LIKE  £XF(-K*C*V/ <BCL12*T)  >  HKERE  V  IS  IK  CK-1 

_  .  AMI  I  JIS  SOME.  1  EMRERAI.LKEx.  LAT_£_R.  HE  IEME£EA1UA£  HfiXEXLES  .. 

ARE  ORE  MULTIPLIED  EV  £k  10  SAVE  ECKE  Cci.FATICNS. 

3ST  =  £CLTZ*TS/ (F*C) 

C  x  SPEED  CF  LIGHT  (CM/SEC) 

H  =  FLANGK'S  CONSTANT  1ERG-SECJ  AFTER  US  1EI1IAL  USE.  KLLTIPLY  ET 
1C*»-Z  TO  GET  XKS  LKITS,  «FICH  Gms  TFt  E.ENUAL  RADIANCES  IK 
_ LERAS...  CF  HATIS . RATHER  THAN  ERGS/ EEC. _ 


A l  r  SORT  (ALOG  (2,  3  )1 
EARTHS  X  EZ£1« 

IS  =  <56.0 
BCLTZ  =  1.JS1E-16 


1  63 

H  x  6.62E2E-2Z 

1  66 

_ 9K.  =  60LTZ/  IK* Cl 

i  65 

8ST  =  TS*ex 

1  66 

_  «  =  n.*i.nt-az 

1  M 

NITR  =  50 

INCLUDING  CLOTES) 


58 


to 

c 

V  10E 

=  V I E  RAT  IONAL  tNERGY  (CM-ll  FF  THE  RAOIATIVE 

61 

C 

TfiANSflTCK 

62 

c 

V I3L 

=  VI  B  RAT  ICKAL  ENERGY  (CM-ll  "F  THE  LONER  ST  A  IE 

63 

n 

v  rao 

s  V  T  Rfl  AT  T  fKfll  OllANTlIM  Ufcfru  fklfB*; 

6<* 

fei 

c 

_ _ £ _ 

THE  EXPRESSION  FOR  THE  PARTITION  FUNCTION 

66 

S3B 

CARO  16  SFECIES  CARD 

67 

at  _ 

;  fiftANT.H  (P  ,  q*  Ci?  ft  \ 

68 

c 

KSOL 

=  ROTATIONAL  CGANTUM  like 

tq 

c 

list 

70 

C 

LST 

=  THE  lcher  STATE  INCICES  (AFGL  like  file  nctationi 

71 

Q 

KTP 

a  i-ItL£S-fcA££  £ P T T  f n y  u  tr.PPH;iTl  =  vnf.Tt  l  =  pnPP^  ct 

72 

c 

VMIN 

=  LOWE*  LIMIT  ON  THE  WAUNIKPE*  CF  ThE  CESIFEO  LINES 

_2J_ 


74 

_zs_ 


76 

77 

78 

_U_ 


8  J 
ni 


82 

_4i_ 


84 

as 


86 

_dZ_ 


88 

_aa_ 


90 

_1L_ 


92 

_3J_ 


94 

_as_ 


96 


98 

_aa_ 


too 

m 


182 

183 


184 

185 


18E 

187 


188 

189 


110 

-XIX. 


112 

iit 


-<m  _=  Ilffia  LIMIT  CK.  THE  NiVENUM-EEA  OF-  Il£_CtSIA£Il.-LltlE.S 


NH-K  RK  ;  4.  ALL  P.  C.  AND  B  PRiNfHFS  ARE  H/AIUATFC 


C 

c 


WHEN  NROl  =  995  ALL  QUANTUM  LIKES  ARE  EVALUATE1" 


DEFAULT  CK  V4 I K  IS  ZERO;  THAT  IS,  THt  L OK EE T -N A VE KLME £ f  LIKES  LN 
TH;  E  1 L E -  A_E L  E VALLLAmi .  HEf  AULI  CN  V«A*  IS  A  NUHEEE  GREATER  T UN 
AAV  CfKCEIvAALE  WA,£NUM8ER  CK  A  LIKE  FILE,  50C"0  CA-1. 


C 

-_C.. 

C 


IAOItf IOUAL  LIKES  can  BE  SELECTED  BY  SETT1KC  BE  TC  f,  C,  OR 
.AKIl.li£QL  IC  THt  OtSIEEQ  LIKE  NUMBER. 


CARD  1C 


IANGEKI  HEIGHT  CARE- --UNIT  1 - LIST -C°^CTE£-  LEAL,  NO  LEF ALIUS 

TANIN  =■  INITIAL  TANGENT  HEIGHT  (»M» 

- -I1H£ - £— E  INAL _ LAHGEHI—EEICEI  mu _ — 


IKTR01  =  examination  interval 


(  *M| 


»1A 


READ  (l,»l  NCL, ISC,  VIBE, VIBL.vIBC 
HRTT-  <4.111  Mfl  .TSn.«TRFrNTRl  tvTRO 


11  FORMAT  <1H1,4X, 'MOL  -  •  ,A3  ,/,5X, 'ISC  =  IECTCFE  CCCE  =',I4,/,5>, 

— -LltlflE.;:  V-lEfiAI  TONAL. EAIERGY  1LH-1J  OF  .  IHE-  IRAK  SI  TICK  t»»£9«3,/»  - 
26X,'VIBL  =  VIBRATIONAL  ENERGY  (CH-1)  CF  THE  LONER  STATE  =  ',Fc.3, 
— 3/.SV»:?8ieC  t-VlBRATiaN-AL  OLANTliH 
4*  =  F  9, 3  I 

_ NBIIFU,..'! _ _ _ __ _ 


CALL  HCLEX  IO_  ESTABLISH 


ASSOCIATED  NI TH  THE  RADIATING  MOLECULE 


CALL  KGLEC (MOL, ISC, MOL NT,DEGV,FRCT , TE  XP1 


WRITE<4,12)MOL>>T,OEGV,PROT,TExF 
12  FORMAT  ABX.'MOLAT  _ 


l'OEGO  =',I3,/,5X, 

Z'E-RCT _ s'.fLui./.SX. 


3'TEXP  =  *  ,F  4,  2 1 


NGT  =  FLOAT 1MOLHTI  /6.02486E+23 


►  to  * . 


iv.- 


121 

122 


123 
:2<< 
125 
1 £i_ 


127 

Jl? 

129 

131 

132 


133 

_L3J, _ 

135 

-1lU _ 

135 

_L1S _ 


139 

19} 

lil 

142 

m3 

19.9 


:«.5 
1 4 1 
:s7 

149 

15.1. 


IF  (9R.E0.  *4  '»  XX( 21  =  XX  11  1 
IFINRGLf EC,999}XX ( 1)  =  XX(1) 
I  =  IMIVK1N1 
N  =  IF 

NSG(1I  =  'VCIGT' 


WRITE (9*211  I,N ,HOL .ISO.UjT ,LST,8R, >X(2  1 ,  FRCL,  *■»(  3  I  ,kSG  (1 1 
21  FGRFtAT  (5X  , 'PROGRAM  Hitt  SE  ARC  F  THE  LINE  FILE  £  r  T  ►  E  E  N  • ,  1 5 .  *  ARC' 
lie,'  CM*  1  FOR  LINES  OF  >,A2,/,5X, 

2'IiCTCFE  =  ',I4,/,5x, 

3'e9NO  s',  IX, 45,*  -  *,1X,45,/,5X, 

4'ERtNCh  -‘.IX. 41.  1X.A5./.5X,  . .  . 

5'LINE  I  =  ',1  2, IX  ,  45, /,  5X, 

5At,'=  LUESH4FE  0°  T  ION  SELECT  EC  ' ) 

WR  I  TF  19,1) 


•1C 


RE4C  ( 1  ,  *  1  T4NIN,  T 4NF,  INTRVL 
WRITE  (L,e 1)  TANIN.TANF .INfRVL 
El  FORMAT  (5X, 'INITIAl.  TANGENT  HEIGHT  (KHS1  *',I£,/,5X, 
T'FINAl  T4FCENT  HEIgHT  (WNS1  =*,I5./,5X, 

i'EXAMINATICK  TFTEA.AL  (KMS1  =*,T5t 


151 

152 


C4RC  2  AICCSPKRE  CARCS - LNIT  2 

TEE  F1KSI  CARP  IS  AN  A  3 -CHAHACIEg  ICEATlfVING  HE  S.STG.E  UUH  IS 


153 

-1S<L- 


155 


SIMPLY  RFJt  4  NC  MRINTcC.  SLCCESSIVE  CAROS  EACH  CONTAIN,  IK  FCR- 
_  M4  T  S  GIVEN  6ELCW  I 


A...  5  AL  I ITLCE  1KPSJ  LI  M£G£P  xAEJLA  ELE  ) 

TRTMP  (41  =  TRANSLATIONAL  TEMPERATURE  (XI 


IE 

F10 ,3 


V  < EOP>  »  ENERGY  V  THE  LC *ER  STATE  (CN-l) 

_ US.  LS  «  U££EJ»  AND  L  OH  FA  S1ATF  IXALSII  TONS,  -RESPECTIVELY 

SCH  ( ERNCh )  •  BMNCN  <P,  0,  CA  R) 

X  *  ISOTOPE  C COE 


W.ITE<<t,91l  UST.IST 

91  ECJLHAI15X  ,  •  AFGX  -LIKE  FILE  Of  1K1E&ESI— E£  (  IRANSIT1CN  ■  , 

1A«,‘  TO 

Mf  I  TEl<t«9-2l _  .  _ 

92  FCRHAT<5X, ‘BRANCH*  ,2X,  ‘LINE  *,  2X ,  *FREC.  CCH-1)*.2X, 

■ma  mu  £ a-eeiii  * ./  x _ 

J!»A*  FOUNTS.  IK£_T0TAL  NUflB-ER  Cf  LINES  TO  eE  EXAMINED 

_ If-AX-i  0.  . . .  . . ..  .  ... 

M*« 

<MI _ _ _ 

:C0  RFA0<3,lU,EN0*lT5IVINIT,z‘l.Z2,Y,u¥,LS,BCH,N,  I 

1111  f  C  RNA.li  F U  .  5»  E 1 J.  X.  F5  .  3,FH..3,JX,AE,3X,AE,14X,H,  13, 4 X,  141 

PROfiRAft  EXAMINES  _ANC  RETURNS  CORRECT  TRANSITION  STATE,  LINETSI 
ANC  BRANCHES)  FOR  THE  GI  YEN  ISCTOFE.  IF  TNE  OCPFlER  CFTICN 
_ tNP.P  =  1)  TS  SFLEOTED.  THE  1  flEr  NT7  LILiKICIHS  ARE..SEI.IQ  ZEJLO  . 

IF<VINJT.IT.VHIN»GC  TO  ISO 
IF I VI  NIT . GT  .  VHAX  >  G  C  TO  1 J 0 
IfdSG.NE.I)  GO  TC  IOC 
IFIUST.NE.LS)  GO  TC  IOC 

IFIBR.NE.  *<*.  AN0.8R.NE.8CH)  GC  TO  100 
IF  1 NRCL.  N£»  N,  AND,  HPCL,  NE.999)CC  TO  190 

JP«*  .  JHFX  *1  . .  ... 

VfJHAX)  =  VINIT 


ALFCJMAX)  =  Z2 

If  1NQF,  EQ.UALF  UN  AX)  =  fl-  fl 

ECP(JNAX)  =  V 

mhchimsaxj  =  ecH 

RCL(JHAX)  =  N 


AT  THIS  FCINT  THE  CATA  HAVE  ALL  BEEN  READ  IN  »nq  STOREC  IN  THE 


ThE  TOTAL  NUMBER  OF  DIFFERENT  LINE-OF-SICH  PATHS  TO  BE  CONSItcREO  IN 

ALL  IS _ (IANF-TANTM/TKTRQL  ♦  1.  THFSF  CASFS  -AJEj.-Q£M£-Itl  EBCI1FS  Of. . 

FOUR  CPOSSIflLY  LESS »  IN  THE  CASE  OF  THE  LAST  £FOlP|  1C  MAKE  TFE 


AF  F  AYS.  THE  LOCP  OC  3UC  IS  EMEREO  ONCE  FOR  EACH  GROLF.  THE  IMEGEP 

.  A  SEAT - HIS  .  SI  CAES  1HC-TANGENT  HEIGHTS  UEJOll  -JitICH  PESLEKEUilZE  -EACH 

PATH.  NREPS  BECOME  S  THE  NUMBER  CF  GPPUFS  MUCH  MUSI  EE  HANCLEC.  IMaV 
.BECOMES  TEEN.UFIIER  CF  CASES -XL  .IhE  CURRENT  JG&OiFU-  . 


NF  =  MCC  IN  REPS  ,  4) 

NREPS  =  ARIES/ A_ 

IF(NR.GT.O)  NREPS  =  NREPS«1 


C»  LCCP  TO  SELECT  GROUPS  CF  IE  OR  LESS!  TANCENT  HFIGHTS.  CALCULATE 

C! _ TtESE.  CASES  SIEUuI  ANilQUSlY  .  KEEN  This  LC  IF  ENCS,  EiECLIION 

C»  TERMINATES. 

Ci _ 

DC  TO.  NF  =  1.F9EPS 


OC  14  G  I  -  1.  N 

_ ttlSELI-  ==  IANIN  ♦  C4»CNR-1)  ♦  I  -  IJ'INTSkL 

IF  (HTSIII  .GT.  TANF)  GO  TC  140 

_ 1KAX  =  I 

140  CONTINUE 


•  SLMRAC  IS  AN  ARRAY  THAT  STORES  THE  SUM  Cl  FACIAnCES  FFCH  HCTAIIONAL 
u! -TRANSITIONS  FCK  A  GWEN  IANGENI  FEIGHI.  F  RESET  IKIS  ARE  AY  IC  ZERO. 

•  ALSO  CALCULATE  THE  ARRAY  OF  PATH  LENGTHS  III  FI-R  each  layer  m 

L* _ AN-CLEACH-FAIH  LIANG  ENT  HEIERIi  IN  THIS  GfiCLP*  ACTUALLY.  ZZ  IS  THE  P-ROtLCT 

:•  OF  THL  TOTAL  NL'MOt  F  DENSITT  AT  THE  ALTITLCE  IF  (JUSTICE  ANC  THE  6E0- 


KM  TO  CM. 

OC  145  I  =  1,1  FAX 
SLitRACLlJi  C.O 
OC  145  N  =  1,  NLYR 


C«»  LCCP  TO  FEEFORN  CALCULATIONS  CN  ALL  SELECTED  ROTATIONAL  L 
C”  QC  COPRECIS  PCS  THE  STINULATEt  EMISSION  A  T  TS  «’«£  K)  . 

C»»  THE  DOPPLER  LINEWICTH  IS _ - _ 

C»»  ALEC  =  SaRT(2»«LOG(2)»€CLT7,TRTMP/(MCLWT*C»C»)«VIMT 


...0C._2ll£..J  =  1.JEAX 

V I  NIT  =  V<JI 

_  ac  =  i.o  -  expT-viNtT/esTi 
OCPP  =  «2»SCRT(2.u»ROLTZ/<eK»kCT»C*C»)«XlMT 


♦  THE  LCCP  CC  150  CALCULATES  EATH-INCEPENOe  KT  QUANTITIES  HriICH 

•  DEPENf  QKLT  ON  AlTITUQE.  ALTITUDES  RUN  FRCH  THE  LCNESI  RECU1SEL  FOR 

•  THE  FIRST  (LOXESTI  LI N £ -OF  - S  1GH T  PATH  TO  THE  HTOh£ST  REQUIRED  FOR  THE 

»  LAST  CNE  IN  THE.  CURRENT  6RCLF .  NHAX  IS  JEL  NLtEEl  CF  SLCM.  RECLIF-EO _ 

»  ALTITUOES. 

* 

•  STSU*  GIVES  THE  LINE  STRENGTH  ENTER  CONDITIONS  CF  LIE  AT  TEMPERATURE 

•  TS.  HUGH  HE  REFER  TO  AS  "STANDARD  conditions**.  the  clantity 

•  QRAT»EXPI21I  GIVES  THE  RATIC  CF  THE  PROBABILITY  CF  FINDING  THE 

•  HC1FCULF  IN  THE  1  QnFR  RAQIAUAG  SIAIE  Al .  AITIIUD£__t  LkE-t-S-S  AflifF _ 

•  AND  TRTHp  CHARACTERIZE  ThE  DISTRIbLTION  OF  POPILATEO  STATES*  T C  THE 

•  PRC8A9KITY  OF  FINCING  IT  IN  THE  L CHER  STAIE  L'NQER  STANDARD  CCNDITIQNS. 

•  THE  RATIC  ll-GANI/CC  CORRECTS  FOR  STIMULATED  PHTSSICN  AT  ALTITUDE 

•  A  ANO  FCR  STANOARC  CONDITIONS.  THE  PRCOLCT  CF  THESE  THREE  FACTORS 

•  GIVES  THE  LINE  STRENGTH  AT  ALT  ITUOE  A,  StAI.  MAMEIM)  GIVES 

THE  CCNVENT  IONAL  A  BSQRPT  TON  COEFFICIENT  A  I_£RICU£NO  LL. _ 

• 

_  NHAX_  =  HTSJIMANI  -  HTST1*  ♦  NLYR 
OC  15U  N  =  1,  NNAX 

.  A_?..MT$U)  ♦  N_  -  1  .  .  . 

GAN  =  exp(-(VINIT-VI8£»/TRTHPIA»N»EXPI-*IEE/VETHF|All 


flagui  *  : 


_Itf  L ££#  tc  lta  calculates,  fcs  ALL  ALT1IL0ES  RE  CL  I* EC,.  THE  VALUES 
OF  THF  VCIGT  OFOFILE  AT  THE  "CISTANCE"  F » C  ►  THF  CENTER  OF  THE  LINE 

INCEXiC  Et  WE  CURRENT  VALLE  Cf  L---IHAI  It..  EVViLI  XF.--1  f  ECU ..  WE _ 

CENTER. 


OC  16  C  N  s  l.NFAX 
A  ;  HIS(  II  ♦  N  -  1 
alfo  =  n«ir(Ai 
X  s  OVVtLi/ALFC 

RAT  r  A2»  »ALF  (  J)»ALC0R  (A)  > /ALFC 

—VGJiAi-—  NNERF  tT.gAM  /ALEC -  -  - - - - 


thf  I  nnp  rr  i;r  I  :  i,IHAA  VAAItS  IKE  TANGENT  HEIGHTS.  THE  PECF1LE  ARGA*, 
N«Z,  IS  also  INITIALIZED.  AnC  the  flag  IS  CHEC*£C  INC  FOSSIBLT  RESET. 


oc  it:  i  =  1. 1 fax 


TFIFLAGII I.E0.1IG0  TO  169 
-A  i -HI  Si  II 

zi  =  mai  »zz<a,d  »vgt  (ai 

_ SIASJ  -THE  Zrlll  IEGRAT10 A  ALCNG  THE  LIFE-OF -SIGH.I  AI  THE  FAR 

HCOIZCF.  TFE  L  CCP  CO  1E5  CARRIES  THE  CALCUIATICN  TC  THE  TANGENT 

-  --HEIGHT  PQ1KT. -IN.  IKXS  LOOP.  ANC  IN  THE  NEXT  CHE,  ZJ  NECCHES  THE 

OPTICAL  OEFTH  CF  THE  CURRENT  SLAB  AT  THE  Fc£CLENCt  IN  QLtSTION. 


nc  it*  n  *  i.nltr 
.A  =  .IUA..I  HTS ill  -  N- 
Z1  *  S<Al»ZZlA,n*iGTU> 


IF  (Zl  «LT.  35.)  Z  Z  *  EXRC-Z1I 


CCHPXEIE  I  Hr.  Z  -  INI  EGBAI ICA  BT  CARR  TING  IT  FRCH  THE  TANGENT  HEIGHT 
POINT  THROLGH  FTLR  SL*eS  TC  THE  OBSERVATION  FCINT. 


nc  17 C  N  «  I.NLTR 


the  LOOP  OC  19 C  RJNS  THROUGH  THE  TANGENT  HEIGHTS  IK  THE  CURRENl 


iijai  mti 


•  •  12  IS  THE  St*  OF  OF  EVEN  TERNS  «THE  ENCS  tf  EACH  PANEL)!  21.  CF 

•  •  IHi  RESULT.  RACUi  ,  MLS  UMJII  CF  NATT/XMIrSIB 

•  • 

lao  dc  m  i  <  uihlx 

THINR1I)  =  0.0 

_ nc  las  M  .  1.K1TR _  _  _ 

A  =  NtTR  «  HTS(I)  -  N 

lES  THINR(I)  T  TUI  KRCI  )  ♦  G3diAJ.*£  (A1*22(A«I) 

THINNED  =  2.0*THIKR(tl  »H*C*VIMT 
RifliXi  «  IfcXNR.iL) 

IF(JT< JI.6E.ltGO  TO  19C 
_ 2?  »  fl.fl _ _ _ 

21  =  L.O 

DC  109_  L  *  2.11.2 

22  *  22  «  NVZCL.II 
Z1  *  il  *  MZlLtX,  II 

1«8  CONTINUE 

22  .  72  .  KW7(11C.  I) _  .  .  . 

RAO ( I I  =  (KVZU.II  ♦  4.9*22  ♦  2.  C* 21  »  N»  2  (10  1 ,1 )  I / 3 .  C 
RAOdl  =  FAC(I>»2.  CE-04*K*C*<IMT 

IF(  (ACS(THINR(I)-RAC(I)  •>  .LI.  (RAC(  Il*.55l  )  JT(JI  *  I 
190  SUNRAOO)  *  SUHRA3  1I)  ♦  RAD  (I) 

•  « 

*  *  END  THE  LAST  LCOP  (INDEX  I)  THAT  VARIES  Th.E  TANC.EKT  HEIGHT _ 


ODD  TERNS 


XF  ONLY  ONE  LIKE  IS  CALCULATED.  PRINT  CUT  IH£  LINESHAfE 


IFtN.P.a.OIN  =  IMAX 

If  UHAX.NE.DGC  TO  195  .  _ 

IF  (K.c0.5)GC  TC  19  4 

«FITEC4.1«liaR.NRQL,iMTStI>  ,1-l.lHAXl  . 

151  FCRHAK///,*  RAOIANCt  (  PHC  T  CN  S/C  H2  *STR  *SE  C*CH  - 1  )  FCR  U 


49S 

_510 _ 

501 


\  MdIHm.'B f  i.ida'i i 


27X,*FPEQ*,2X,4(2X,I9,«  KH')> 
192  FCRHAT  (FX1.4,  2X.4E14.5) 

HR  ITE  (4 . 2) 

DO  193  L  =  1.1C1  . . 

IF<NVZ<L,M  .EC.0.  0  )GO  TO  194 
T 


193  HFITE(4.  192)2  l.(NV2(L.  II  .  1  =  1 .  N) 

ISA  NFITE  <•..!» _ 

HRITE  (4.140  KHTSfl  ).HTS  (I)  .  I=1»IHAX) 
HP;TE<4i2I  . _ 


195  HFITE(4,19E)BRKCHLJUR£LX4IillfcXfcBXXJ  .1=1.  IMAX1 
IF(N.LT.5I»RITE(4.  1971  (RAO(I)  .I=1.N) 


SI* 
A1A 
sit 
All. 
sit 
A1A 
520 
ALL 
522 
A2A 
52. 
A2A 
52t 
522 
52  ( 
A2A 

531 
AIL 

532 
AAA 
S3* 
AAA 
33t 
AAL 
531 
AAA 
5*0 


lit  FC9NAT(5X,M.I*,.(2X,E12.5.12»>> 

- 11?  FTIHIT  MH9.9W.M  UI.£A2.»5I  I  -- . 

imjt(  j)  .ei.oun  Jl  s  ; 

XAA _ _ _ _ 

210  CONTINUE 

-LAS _ _ _ _ _ 

C»*  END  TNE  LOOP  ( INOE »  J»  THXT  CFOOScS  OIFFfRENT  (INES 
XA * . . 

C*  . . . . . . . 

xa _ _ ....  ....  .  . . .  _  _  _ 

MR ITE (*» 1) 

- 1/  mmuAUu  NtriF  umXaiajjuuu  mmti  tit  tm,  !►.»■ _ 

201  FCRNAK*  TOT4L  FCR',1*.*  LINES  =  *,.(  E12  .  S,  l.v>» 

300  CONTINUE  *  .  . . 

C*  END  TFE  LOOP  ( INOE  2  Nil  SELECTING  CIFFERE NT  GROUPS  OF  IINc-CF-SIGHT  p«1hS 

XA _ _ _ 

C 

X _ _ _ _ _  _  _  ..  _ _ 

99»  IFUfMX.tC.OIMRm  (*,  9901 

- 990  XC9NAT I///, «  .....  an  itnfc  m.>n  .«»»«»,./#/>  _ _  _ _ 

CLOSE  UJMT«1I 

- CIOSF  (IJKTI.?I _ 

CLOSE  (UNI I *3  I 

- Cinsf  (UMI.NI _ _ _ 

ENC 


\  »  ■>.*. .''MiiV.  ■  ■I'.Vw. •- ■'  '•  *■  •-  -  » i.  •  •  ...  mLt  t. 


V.SV 


THE  Ht  THCC  IS  cue  TO  RT8ICKI.  THE  FGRTRA I  LISTING  NAS 
“FORNATICN  Of  LINE  ANO  CONTINUOUS  SPECTRA,”  SPFCIAl  EFFORT  MI, 


kin  if  I  w  Wi  I  ill  Ja  Ml 


C  THE  HETHCC  OF  _  CAUtULAllQN  IS  EESCRifiEE  It  fl.  H.  AFftS-LEORG 
C  ANO  R.  «.  SICHCLLS  ■  “EPISSICN,  ABSORPTION  ANO  TRANSFE  R  OF 

_C _ RADIATION  IN  HEATEC  ATMOSPHERES,”  PERGAHCN  PRESS. KEN  TORN. 

C  119721  -  PP.  235-Z37. 

r  r  .pi  n  7 


CINENS1CN  Cll.'ll 

. c*ta  ntry.kmax.rh.pi.czjii:  /t,is,!.Ci:.:u$i2es!!(;t 
1  5.6<.1A9**354776E-01,  A.9753*E10E258JE-02/ 

IF  IMRT  ,EO.  1)  GO  TO  170 
11C  CONTINUE 

_ 3  f  £ILK2»Ut _  _ 

IF  <«  .EC.  O.tl  GO  TO  160 
_  .  IF  ilx.  ♦  t|  ,CT.  25.81  GC  TC  190 
nc  CONFUTATION  FOR  GENERAL  CASE 
_  _  A1  .*  RH»J 

A2  =  A»A 

if  l A  .lt.  a. it  G£  IC-128 

Z  =  C£XP«CW>LX(-01»A,C1*XI) 

VR-RF  =  C.O 
GO  TO  1JC 
120  CONTINUE 

Z  =  CCCS  UNFLX((U»X,01»AI  ) 

_ NWEP-f  S  C2«£Xf(A2  -  K«  XI  «CCS  (2.  £?A».X1 . . . 

130  CONTINUE 

_ £1  r _ <  1_»£  v.RiAUZ  >»•<*<).  5  •**- 

ez  =  -ainagcz) 

_ _ S  *  -0.5*(FLCATINfl)  ♦  RH* x) 

T  =  S»S  4  0. 25*RH  SO*A  2 


T  =  T  4  I  4  0.25 

-iii-'-W  - 

ei  =  ai  -  at 
£2  5  -B2 

IF  IT  «GT.  2.5E-12I  GO  TC  140 


I’M; 


GO  TO  15  C 
INC  XCNTUNUE 

NHERF  =  8NERF  4  CIN»*(82*S  4  9D/1 
150  CONTINUE 
155  CONTINUE 


«  =  -NP1 

QC  lh_  Hal.  HPUPl 
K.  -  K  ♦  1 

CtUl  i  CJ^EXRi-f  1.CAIXK.*JCl/fcl--SC) 
let  CONTINUE 

*tln:  *  sort  ulog  <2. 3 )  > 

CO  TC  lie 
1=0  CONTINUE 


USE  4S:rNEICTIC  EXPSNSICN  PC*  CCNPI.C.X  ERROR  FCh'-TICT 

fog  LARGE  »  Ut  1. _  .  - _  _ 


-JUl  S  -2.E 
CENP  =  fl 
C£M  =  -J 
OC  20.  I  *  1»  E 

n.EN  -  flh/trFNS»ilF  NR  t-  rFM  »r,  EMI  I  - 
C;np  =  CEX'OEXR  »  s 
— C£HI  _=  -CENT  *£EN-  -  X 
4N  =  SN  -  Q.S 


TEN  =?  02/1  C£NR*OE  XR  ♦  OtNI'CEM) 

uyCoc  «  o t  ■  iioioik  r  jmi  r  Cud 


fi  tTUFN 


SU8R0UTIKE  HOLeC(NCL,  I SO,PCL«T,CEG X, PROT , TE XP ) 


THIS  SURR  CLTI RE.  GIXtN  THE  HCLECLL  £  COOE  (t>GL)  AKO  THE  ISOTOPE  COOE 


HCLECtlLE  IK  QUEST ICNI  NOLkT,  THE  MOLECL'l AR  kETGHT 

n:r l  tu  c  rcrr a  p  ■  -  - - 


PROT.  THE  EXFONERT  OF  ’t^FEMURt  IK  THF 

-  -  -  -  hCIAIICNAL  EAB.I1UXN  FLKCTION . . 

TEXF,  THt  EXFCKEM  OF  TFRFERATLRE  IK 


IX  l*>  TO  *  1 0  K I  IF  A  SINGLE  LE  vE  L  AT  1J?C  CM-I  NITF  A  CE- 
generact  of.  .5.  is  -Chose .  _ _ _ 


QEGX  1S-Z-  .fO a.  LLHEAH  .-t&LA icmo  PQLfCULES..  1  -0 TttEJUXSX -LXC-OU. L 
PKCT  IS  1  FOR  ALL  LINEAR  HCL£CUL=.S.  1.5  CTHERNISE 


INTEGER  DECS 

CHARACTER**  MOLCOO  (T)  ,KOL 
- OmU&HU.  Jt  CL  MULE  UHQLISQIH 

OATA  CH0LCCCCII»I=1»X|/*M2C*.  *C02*,*03  *,*K20*,*C0  •,*CH.»*,,02  V 

_ Oil  A.  1-MJ1  I?C1I  T  UUUIUm  ■  F  :  i ,  UK  -j.  _ _ _ _ _ 

OATA  (MCLRG  T<  11,1=  1,7)  /  la,  1, A ,  L8,  AN  , 2  6,  16, 11/ 

NCLNT  =  o 

-  -  IC  *  i  .  . 

PR  QT  =1.5 


TEXP  =  0.5 

_ Of.  10.  1_-S_  UI- 

ISC  IF  (POICCCUI.  EC.NOLIIC  =  I 

_ LE  irr-FQ-niBgruoM  .. 

IF(IC.tQ.l.CR.IC.£Q.3)G0  TO  12 C 


IF(IC.NE.?»PRCT  =  1.0 

IF  (IC.EQ.  7IOEGV  =  5 

jigitF-  =.q<L_a  mx scoot _ 

IMNOIFFIIJO,  1A0.150 


H2  =  INDIFF  -  100»M1I/10 


ISO  *  ISCTOPE  COCE  >  62% 


2K.3SC  198.393 
_2aa.39n  „x9(..o69 
236.636  191. <<66 

-21U.625...  X16.336 
202.520  196.632 


19«. 610  127.799 
X96.65S  12  Hw.hii 
191.710  171.720 


190.780  165.253 


.276623£-C<. 
•23SS2J6-U.. 
.2OC670E-O4 
.12 162  a  E  --06 
.157219t-J« 


.1028735-3* 

.8760906 -05 
•  7  6  029  3  S-.5 
•622aXj6-M»5 
•52510CE-35 


lae.aio  160.017 

t«e.»7C  156.397 

1*6.870  156  .  <.0* 


187.360  155.926 

188.330  156.292 


.16  6500  E-35 

,2  56  310  E-35 
_^1UUI-15 
.1611905-05 


•  12  7190  £-35 
_  -  -10  6633  6-35 

•  89I753E-U6 


.303800E.12 
.  2612(067 1 2 
. 22O730E* 12 
.192  7-506*12 
.  165  32  0E* 12 


. 123913E* 12 
.  10-3C5CE*  12 
•  8  79  Cl  QE.  11 
.  7471406*11 
. 635C50E* 11 


.*  55820t*ll 

. 319363E«11 
.26ZX.CX*  IX 
. 223800E* 11 


. 156530E* 11 
..  .  130  6-5H6.lt 
. 10S380E* 11 


195.080 


200. 2C0 


205. 31C 
20.6.41? 
213. 52C 


223. 98C 


.9157006-07 


VfM.rin 


*  9C38  706*  C9 


141^^4^X111:1 


113 

240.000 

HC.  *16 

.7011306-0  7 

.  54  72706*  09 

111 

_ 2-51.1*0  - 

X\kxZ51 

.015970-6-0  7 

.4275506*05 

»* 

112 

264.000 

187.580 

.54621*6-07 

•  334*906*  09 

.* 

113 

_ 276.1 1£_ 

19 C. 715 

.49C53QE-07 

•  26b  7405*09 

114 

280.030 

193, 65*. 

.43892  JE-C7 

.  2125(06*  09 

M'lnTH 


-U2L _ 125.  (Lit 


336. OCC 

_ 398.00  0 

360.000 

_ 371.995 

382. 89C 
395.335 
9P5.78C 

_ Alt./lO 

427.69  0 

_ 5lti.t»5 

448.950 
**56.8  6  0 
469.270 
478.305 
487.340 

_ *31.1-7-5 

505.410 


Hi _ 514x5  is 


523.490 
_5  32,  525 
541.  56  0 
551.595 
559.638 
ife7.ux_ 
574.560 
-581x055 
589.530 


604.490 


152 

153  _ 

159 


634. 390 


646.770 
652. 96C 
655.150 


200.369 
203.701 
206.862 
209.852 
212.375 
_214x-411.- 
216 .214 
217.124 
215 . 243 
2  2  Gx_5  7  9 
221.845 
222.954 
223  .  9"c 
224. 727 
225.533 
226.346 
227**53 
-117x442- 
228.177 
-224x6*1 

225.068 
225. *27 
225.723 
■225.955 

230.139 

-Zll-llk 

230.368 
-230 .423 
230.441 


r»»itiT 


'  (0 . 385 
23  C.J14 
230.230 
.231.1X9- 
225.979 


220.636 

-219.995 

225.242 


iin  c  «i 


.3583702-U7 
.32652*  6**7 
.297530c -0  7 
.2729976-07 
.2505005-07 
-.21127  0  I-*  7 
,21 352  0  £ -0  7 
.1981605-07 
.16391*6-07 
, 1714505**7 
.15983  0  5-07 
,19S5/.li-J7 
.1  359*06-07 
.1314336-1,1 
.123410  E-u7 
.11(2306-07 
.1054706**7 

.10217.04 -.07 
.9761405-08 
.92397-36-04 
.87459*6-08 
.8256  20  6-09 
.78696*  6-0  8 

.7108706-38 
•  67192*6-0  8 
.644400  *-68 
.619930L-U9 
. 5860 9J 5-09 


WWiU 


.5345*06-08 
.5112106-08 
.*889305-08 
.968150  6-0  9 
.4482806-08 


TTI'l.TIlin'M 


.41367*6-04 

”.3817306-08 


.  1  40  76  06*  C9 
.  1162206*  09 
.5595506*08 
.8068706*08 
.  6717(05*  08 
.5612502* 08 
.4  82  5406*08 
.  412550E*  08 
.  3534306*  08 
.  305  38  06*  08 
.263(806*08 
. 229S4i)£+  Cl 
.20*  3606*  C8 
.  17  (  7006*  38 
.  1540  7*6*  08 
.  1*62306*  08 
.  U3  4  5  06*  08 
,lH52fi£+M 
.9532(06*07 
,8»H306*C7 
.  75607*6*  C 7 
.  6  7706  06*  07 
.6063006*07 

-  54*5. of + o y 

.  4  85  75  0t*  C7 
. 442(306*  07 
.4003806+  07 
. 3629606*27 
.  329*306*  07 


IQifi 


.  2717506*  07 
.  2475606*  07 
.2255206**7 
•  2  0  5  510  6*  07 
.  1  880006*  07 


m  wn:r4'l  xzi  wm 


.  15905  06*  07 
_ul912906»C7_  . 
.  1345606*  07 


40N 


716. 8C2 


727.058 


737.315 


747.570 


226 . J55 


225.356 


224. 662 


223.955 


223.  199 


.259760  6-08 


.2  41730  E-l*  8 


•22495u£-08 


.2063305-08 


.1960206-08 


-u8 


•  6  02  740  E*  06 


.  5  184206806 


.4458906806 


.  3  8  3510E8  C6 


.  3339106*  06 


.  290  72  OE*  06 


176 

771.948 

221. 757 

.1718846-08 

. 253120t*C6 

K'. 

_ 122 _ 

_ 275*011- 

221^382 

.161333  6-08 

.2361901* 06 

7»'J\ 

in 
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